HALFTONING METHOD AND APPARATUS, AND 
COMPUTER-READABLE RECORDING MEDIUM IN WHICH 
HALFTONING PROGRAM IS RECORDED 
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BACKGROUND OF THE INVENTION 

1. Field of the Invention: 

The present invention relates to a halftoning 
10 method and apparatus for converting a multilevel input 
image, such as a continuous-tone color image, into a 
binary image, andalsotoa computer- readable recording 
medium in which a halftoning program is recorded. 

15 2. Description of the Related Art: 

For reproducing a tone-graduation image (e.g., 
multilevel image and continuous-tone color image) on 
a printer that is dedicated to binary outputting, it 
is necessary to convert values of the multilevel image 

20 into binary values . The binary-type printer reali zes 
reproduction of tone graduation by modification in 
terms of dot area or dot density. In the following 
description, the terms "binarizing" and "halftoning" 
are used each as a word having a meaning "converting 

25 a multilevel image into a binary image". 

Consequently, in such halftoning, it has been 
required to make moire and artifacts, which does not 
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appear in the original image (multilevel image), 
unconscious and also to reproduce the tone graduation 
of the original image with the best possible fidelity. 

The conventional halftoning technology is 
exemplified by the following: 

(1) Japanese Patent Laid-open Publication No. 
HEI 6-070144 (Application No . HE I 5-142746) discloses 
an attempt to realize a moire- or artifact-free binary 
image by carrying out the error diffusion method along 
"a random spatial filling two-dimensional curve" 
throughout the entire original image. 

(2) Japanese Patent Laid-Open Publication No. 
HEI 6-006586 discloses an alternative attempt to 
realize a moire- or artifact-free binary image by 
minimizing energy functions that are calculated from 
values of the original image (multilevel image) and 
those of a binary image, utilizing the "simulated 
annealing method". As this method two technique are 
currently known; one technique is to directly 
converting values of the original image, and the other 
technique is to previ ously obtain 256— grade dither 
patterns prior to converting values of the original 
image . 

However, these conventional technologies would 
encounter the following problems: 

(1) As to the technology disclosed in Japanese 
Patent Laid-Open Publication No . HEI 6-070144 , because 



an entire image is processed by the same algorithm, 
moire or artifacts would be inevitable with the 
converted binary image. In a patterned image such as 
a checkered one, in particular, horizontal bands tend 
5 to appear significantly apparently after halftoning. 
Under the method of "the random spatial filling 
two-dimensional curve", since scanning is performed, 
during conversion, in a single continuous stroke-like 
manner, it is impossible to take parallel calculation, 
10 thereby requiring a huge processing time. 

(2) As to the technology disclosed in Japanese 
Patent La id-Open Publication No. HE I 6-006586, because 
huge time is needed in calculation during conversion, 
this technique would be less useful. As another 

15 disadvantage, in reduction of a very gentle 

tone-gradation image, a peculiar artifact would occur. 

(3) In the above-mentioned conventional 
technologies, the multilevel pixel value of an object 
pixel (noteworthy pixel) to be halftoned is converted 

20 into a binary value, and an error accompanying the 
conversion is diffused to other pixels adjacent to the 
noteworthy pixel. So, in such technologies, if the 
pixel value of the noteworthy pixel is outstanding due 
to an image defect, a significant error risen from the 

25 outstanding value is also diffused to unprocessed 
adjacent pixels. For example, the moire that is 
latently contained but not visually apparent in the 
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original image, which moire is due to periodic 
structure typical of an input device (scanner, digital 
camera, etc) , could be enhanced by the error diffusing 
so that it becomes visually apparent after halftoned. 

5 

SUMMARY OF THE INVENTION 



With the foregoing problems in view, one object 
of the present invention is to provide a halftoning 
10 method which prevents the moire that is caused by sudden 
shifts between pixel values due to image defects in 
the original multilevel image and which also reduces 
moire and other artifacts surely to a minimum by a simple 
technique . 

15 Another object of the invention is to provide 

an apparatus for carrying out the above-mentioned 
halftoning method. 

Still another object of the invention is to 
provide a computer-readable recording medium in which 

20 a halftoning program is stored. 

In order to accomplish the above first-named 
object, according to a first generic feature of the 
present invention, there is provided a halftoning 
method of converting a multilevel input image into a 

25 binary image, comprising the steps of : (a) calculating 
the multilevel value of a given noteworthy pixel of 
the multilevel input image, as an estimated value of 
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the noteworthy pixel, based on the multilevel values 
of pixels other than the noteworthy pixel; and (b) 
converting the estimated multilevel value of the 
noteworthy pixel into a binary value. 
5 The above second-named object is accomplished 

by a halftoning apparatus for converting a multilevel 
input image into a binary image, comprising: an 
estimating section for calculating the multilevel 
value of a given noteworthy pixel of the multilevel 

10 input image, as an estimated value of the noteworthy 
pixel, based on the multilevel values of pixels other 
than the noteworthy pixel; and a binarizing section 
for converting the estimated multilevel value of the 
noteworthy pixel into a binary value. 

15 The above third-named object is accomplished by 

a computer-readable recording medium in which a 
halftoning program for instructing a computer to 
execute a function of converting a multilevel input 
image into a binary image is recorded, wherein the 

20 halftoning program instructs the computer to function 
as the following: an estimating section for calculating 
the multilevel value of a given noteworthy pixel of 
the multilevel input image, as an estimated value of 
the noteworthy pixel, based on the multilevel values 

25 of pixels other than the noteworthy pixel; and a 
binarizing section for converting the estimated 
multilevel value of the noteworthy pixel into a binary 
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value . 

With this construction, at the conversion of the 
multilevel value of the noteworthy pixel into a binary 
value, the value to be converted is not the value of 
5 the noteworthy pixel itself but the estimated value 
of the noteworthy pixel, which was calculated based 
on the values of pixels other than the noteworthy pixel. 
Therefore, even if some pixels have outstanding pixel 
values due to image defects, significant errors risen 

10 from the outstanding value can be prevented from being 
directly diffused to non-binar ized adjacent pixels. 

The halftoning method may further comprise the 
step (c) of diffusing a possible error, which has 
occurred in binary value with respect to the noteworthy 

15 pixel, to multilevel pixels adj acent to the noteworthy 
pixel by a technique. A possible error, which has 
occurred in binary value with respect to the noteworthy 
pixel, may be diffused to the pixels used to calculate 
the estimated value of the noteworthy pixel in step 

20 (a) . 

According to the halftoning method and apparatus , 
and the computer-readable recording medium of the 
present invention, it is possible to guarantee the 
following advantageous results: 
25 (1) at conversion of the multilevel value of the 

noteworthy pixel into a binary value, the value to be 
converted is not the value of the noteworthy pixel 
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itself but the estimated value of the noteworthy pixel, 
which was calculated based on the values of pixels other 
than the not eworthy pixel . Accordingly, even if some 
pixels have outstanding pixel values due to image 
defects, significant errors risen fromthe outstanding 
value are prevented from being directly diffused to 
non-binarized adjacent pixels. It is thus possible 
to minimize moire and other artifacts that did not 
appear in the original (multilevel) image. 

(2) A two-dimensional digital filter 
facilitates the calculation of the estimated value of 
the noteworthy pixel. 

(3) With a two-dimensional digital filter 
dedicated to enhancing the profile, such as a Laplacian 
filter or a Prewitt filter, in calculating the 
estimated value, it is possible to enhance profiles 
in a halftoned image, thus forming a well-defined 
image . 

(4) Since an error diffusion technique is 
changed to another as the scanning of the pixels of 
the multilevel input image progresses, it is possible 
to minimize moire and artifacts in the binary image, 
which did not appear in the multilevel image. 

(5) Since the error diffusion technique is 
changed from one to another upon detection that the 
noteworthy pixel is a pixel constituting part of the 
profile of the multilevel input image, it is possible 



to localize the effects (i.e., moire or other 
artifacts) that have been caused by the change of the 
error diffusion technique to dissolve into the profile, 
thus minimizing moire and artifacts in the binary image, 
which did not appear in the original (multilevel) 
image . 

(6) The change of the error diffusion technique 
for every pixel of the multilevel input image would 
also minimize moire and artifacts in the binary image, 
which did not appear in the original (multilevel) 
image . 

(7) Upon judgment that the noteworthy pixel is 
a pixel constituting part of the profile of the 
multilevel input image, values according to the 
possible error, which has occurred in the binary value 
with respect to the noteworthy pixel, are added to the 
values of the unscanned pixels along the profile in 
the detected direction as an exceptional process. 
This would also make it possible to minimize moire and 
artifacts in the binary image, which did not appear 
in the original (multilevel) image. 

(8) Since the profile value of the noteworthy 
pixel is calculated based on both the multilevel value 
of the noteworthy pixel and those of the adjacent pixels 
and then compared with a predetermined value, it is 
possible to facilitate discriminating the profile. 
Further, the profile discrimination can be carried out 



simultaneously with the scanning of the multilevel 
image, thus speeding up the profile detection. 

(9) With a two-dimensional digital filter 
dedicated to enhancing the profile, such as a Laplacian 
filter or a Prewitt filter, in calculating the profile 
value, it is possible to facilitate discriminating the 
profile . 

(10) The profile value can be directly 
calculated by making addition and subtraction 
individually on the multilevel values of the noteworthy 
pixel and the adj acent pixels . Sincetheprofilevalue 
can thus be calculated without making any 
multiplication or division, it is possible to decrease 
the load of the calculation, thereby speeding up the 
calculation of the profile value. 

(11) Since the error diffusion technique can be 
changed to another that is selected in a predetermined 
order from various different error diffusion 
techniques, it is possible to minimize moire and 
artifacts in the binary image, which did not appear 
in the original (multilevel) image. 

(12) Since the error diffusion technique can be 
changed to another that is selected at random from 
various different error diffusion techniques, it is 
possible to minimize moire and artifacts in the binary 
image, which did not appear in the original 
(multilevel) image. 
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(13) Partly since the occurred error is 
proportionally distributed to the unscanned pixels 
adjacent to the noteworthy pixel in accordance with 
the predetermined weighting pattern, and partly since 

5 the error diffusion technique can be changed by 
changing the predetermined weighting pattern to 
another, it is possible to change the error diffusion 
technique easily and stably. 

(14) If two or more multilevel input images to 
10 be halftoned have a substantially identical profile, 

a discrimination is made, for only one of the plural 
multilevel input images, whether or not the noteworthy 
pixel is a pixel constituting part of the profile. And 
the result of the discriminating can be used in 

15 halftoning the remaining multilevel input images . It 
is therefore possible to halftone the plural multilevel 
input images without making this discrimination for 
the remaining multilevel input images, thus improving 
the practical speed of halftoning. 

20 Other objects and further features of the present 

invention will be apparent from the following detailed 
description when read in conjunction with the 
accompanying drawings. 

25 BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram showing a functional 
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construction of a halftoning apparatus according to 
one embodiment of the present invention; 

FIG. 2isablockdiagramshowinga computer system 
on which the halftoning apparatus of FIG. 1 is realized; 
5 FIGS. 3A and 3B are diagrams illustrating the 

manner in which the pixel value of the noteworthy pixel 
is estimated; 

FIG. 4 is a diagram illustrating the manner in 
which a two-dimensional digital filter dedicated to 
10 enhancing the profile is used; 

FIGS. 5A and 5B are diagrams each illustrating 
a Laplacian filter; 

FIGS. 6A and 6B are diagrams each illustrating 
a Prewitt filter; 
15 FIG. 7A is a diagram showing non-binari zed pixels 

neighboring the noteworthy pixel . FIG. 7B is a diagram 
showing a two-dimensional digital filter applied to 
the pixels illustrated in FIG. 7A; 

FIG. 8Ais a diagram showing non-binarized pixels 
20 adj acent to the noteworthy pixel . FIG . 8B is a diagram 
illustrating a two-dimensional digital filter applied 
to the pixels illustrated in FIG. 8A; 

FIG. 9 is a diagram illustrating directions in 
which profiles extend; 
25 FIGS. 10A through 10D are diagrams each showing 

a weighting pattern according to an error diffusion 
technique ; 
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FIG. 11A is a diagram showing a weighting pattern 
according to Javis, Judice, and Ninke. FIG. 11B is 
a diagram showing a weighting pattern according to 
Stucki ; 

5 FIG. 12 is a flowchart illustrating the operation 

of the halftoning apparatus of the present embodiment 
to change the error diffusion technique (weighting 
pattern) across the profile; 

FIG. 13 is a flowchart illustrating the operation 
10 of the halftoning apparatus of the present embodiment 
to change the error diffusion technique (weighting 
pattern) for every pixel; 

FIG. 14 is a flowchart illustrating the profile 
detection and the binarizing process carried out in 
15 the halftoning apparatus of the present embodiment. 

FIG. 15 is a flowchart illustrating the operation 
of the halftoning apparatus of the present embodiment; 

FIG. 16 is a graph showing the relation between 
the wavelengths and the luminance for magenta; 
20 FIG. 17 is a graph showing the relation between 

the wavelengths and the luminance for green; 

FIG. 18 is a graph showing the relation between 
the wavelengths and the luminance for blue; and 

FIG. 19 is a graph showing the relation between 
25 the wavelengths and the luminance for composite black 
(equal in RGB levels) . 
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DESCRIPTION OF THE PREFERRED EMBODIMENT ( S ) 



One preferred embodiment of the present invention 
will now be described with reference made to the 
5 accompanying drawings. 

FIG. 2 is ablock diagram showing a computer system 
which realizes a halftoning apparatus according to one 
embodiment of the present invention. As shown in FIG. 
2, the computer system of the present embodiment 
10 includes a CPU 10, and a bus line 12 connected to the 
CPU 10. 

A ROM 14 and a RAM 16 are connected to the bus 
line 12. Also a keyboard 20, a mouse 22, a monitor 
(such as CRT, LCD, PDP) 24, and a magnetic disc 26 are 

15 connected to the bus line 12 via an input/output 
interface (I/O) 18. Still also a printer 30 is 
connected to the bus line 12 via another input/output 
interface (I/O) 28. 

And an application program is stored in RAM 16 

20 which realizes a binarizing section 51, an error 

diffusing section 52, an error diffusion technique 
changing section 53, a pixel-on-prof ile detection 
section 54, a direct ion-of -prof ile detection section 
55, and an estimating section 56 on the CPU. As it 

25 executes the application program read out from the RAM 
16 via the bus line 12, the CPU 10 functions (described 
in detail later ) as the binarizing section 51, theerror 
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diffusing section 52, the error diffusion technique 
changing section 53, the pixel-on-prof ile detection 
section 54, the direction-of -prof ile detection 
section 55, and the estimating section 56, thus 
5 realizing the halftoning apparatus of the present 
embodiment . 

The program, which instructs the CPU 10 to 
function as the binarizing section 51, the error 
diffusing section 52, the error diffusion technique 

10 changing section 53, the pixel-on-prof ile detection 
section 54, the direction-of -prof ile detection 
section 55, and the estimating section 56, is 
previously recorded in a computer-readable recording 
medium in the form of a floppy disc or a CD-ROM. In 

15 use, the computer reads out the program from the 

recording medium, and stores the read-out program into 
an internal storage device or an external storage 
device. Alternatively the program may be stored in 
a storage device, such as a magnetic disc, an optical 

20 disc, or a magneto-optical disc, whereupon the program 
may be provided from the storage device to the computer 
over a communications channel or path. 

When realizing the functions of the binarizing 
section 51, the error diffusing section 52, the error 

25 diffusion technique changing section 53, the 
pixel-on-prof ile detection section 54, the 
direct ion-of-prof ile detection section 55, and the 
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estimating section 56, the program newly stored in the 
internal storage device (RAM 16 in the present 
embodiment) is executed by a microprocessor (CPU 10 
in the present embodiment) of the computer. 
5 Alternatively the computer may directly read out the 
program from the recording medium for this execution. 

In the present embodiment, the term "computer" 
is hardware and an operation system inclusive, meaning 
hardware operable under the control of an operation 
10 system. Assuming that the application program solely 
operates hardware in the absence of an operation system, 
such hardware itself is regarded as an equivalent to 
a computer. Accordingly hardware includes at least 
a microprocessor, such as CPU, and a read-out device 
15 for reading out a computer program recorded in a 
recording medium. 

The application program includes program codes 
for instructing the computer to function as the 
binarizing section 51, the error diffusing section 52, 
20 the error diffusion technique changing section 53, the 
pixel-on-prof ile detection section 54, the 
direction-of -prof ile detection section 55, and the 
estimating section 56. Alternatively, part of these 
functions may be realized by the operation system 
25 rather than the application program. 

Further, the recording medium of the present 
embodiment may be in the form of one selected from 
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various kinds of computer-readable media that are 
exemplified by a floppy disc, CD-ROM, a magnetic disc, 
an optical disc, a magneto-optical disc, an IC card, 
a ROM cartridge, a magnetic tape, a punch card, an 
5 internal or external storage device (memory such as 
RAM or ROM) of a computer, and printed matter printed 
or otherwise labeled with codes such as bar codes. 

The computer system of the present embodiment 
shown in FIG. 2 functions as the halftoning apparatus 

10 which converts a multilevel input image into a binary 
image. As shown in FIG. 1, the computer system 
functions as the binarizing section 51, the error 
diffusing section 52, the error diffusion technique 
changing section 53, the pixel-on-prof ile detection 

15 section 54, the direct ion-of -prof ile detection 
section 55, and the estimating section 56. 

FIG. 1 is a block diagram showing a functional 
construction of the halftoning apparatus of the present 
embodiment. FIGS. 3A and 3B are diagrams each 

20 illustrating the manner in which the pixel value of 
the noteworthy pixel is estimated. 

The halftoning apparatus of the present 
embodiment first expands image data of a multilevel 
input image on the memory (RAM 16) or the like. The 

25 estimating section 56 calculates the value 

(hereinafter also called "estimated value") of the 
noteworthy pixel based on the multilevel values of 
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pixels other than the noteworthy pixel, and the 
binarizing section 51 then binarizes the noteworthy 
pixel based on the estimated pixel value. 

In the present embodiment, each pixel to be 
5 binarized is defined as having a pixel value in a range 
of from 0 to 255, having a pixel value of either 0 or 
255 after halftoned (binarized) . This definition 
would never affect the generality of description. 

The estimating section 56 calculates the pixel 

10 value of the noteworthy pixel (reference character A 
of FIGS. 3A and 3B) based on the multilevel values of 
pixels other than the noteworthy pixel A. 
Specifically, a two-dimensional digital filter is 
applied to non-binar ized (multilevel) pixels, other 

15 than the noteworthy pixel , included in a predetermined 
area. One of the following three types of techniques 
(1) through (3) is selectively used to calculate the 
estimated value of the noteworthy pixel A. 
(1) First Estimating Method: 

20 FIG. 3A shows one technique of estimating the 

pixel value of the noteworthy pixel A in which a common 
two-dimensional digital filter is used. The pixel 
value of the noteworthy pixel A is estimated based on 
the values of non-binar i zed pixels in a predetermined 

25 area (2x2 pixels in FIG. 3A) that centers around a 
position predetermined distance (2><6 pixels in FIG. 
3A) apart from the noteworthy pixel A. 
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Using a matrix 100a (for example, 2><2-pixel area 
in FIG. 3A) as of FIG. 3A, a two-dimensional filter 
of FIG. 3A scans 2x2 pixels in the predetermined area 
while making calculations with the pixel values in the 
5 matrix 100a. Specifically, after multiplying each 
value of the pixels in the matrix 100a by 1/4, the 
results of the multiplications are added to obtain the 
pixel value of the noteworthy pixel A as the estimated 
pixel value. 

10 (2) Second Estimating Method: 

FIG. 3B shows another technique of estimating 
the pixel value of the noteworthy pixel A in which a 
digital filter dedicated to enhancing the profile is 
used. The pixel value of the noteworthy pixel A is 

15 estimated based on the values of non-binar i zed pixels 
in a 3*3-pixel area (predetermined area) that centers 
around a pixel located 5.5 x 6 pixels (predetermined 
distance) apart from the noteworthy pixel A. 

Using a matrix 100b (for example, 3><3-pixel area 

20 in FIG. 3B) as of FIG. 3B, a two-dimensional filter 
of FIG. 3B scans 3><3 pixels in the predetermined area 
while making calculations with the pixel values in the 
matrix 100b. Specifically, the digital filter 
dedicated to enhancing the profile is applied to each 

25 of the pixels in the predetermined area to estimate 
the pixel value of the noteworthy pixel A. 

In use of the digital filter dedicated to 
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enhancing the profile, a Lapracian filter (FIGS. 5A 
and 5B) or a Prewitt filter (FIGS. 6A and 6B) is 
preferred. FIG. 3B shows an example in which the 
Lapracial filter is employed. 
5 A manner will now be described in which the pixel 

value of the noteworthy pixel A is estimated using the 
digital filter dedicated to enhancing the profile. 

FIG. 4 illustrates the manner in which a 
two-dimensional digital filter dedicated to enhancing 
10 the profile is used. FIGS. 5A and 5B and FIGS. 6A and 
6B respectively show two-dimensional digital filters 
dedicated to enhancing the profile: FIGS. 5A and 5B 
illustrate Laplacian filters ; FIGS. 6Aand6B, Prewitt 
filters . 

15 When the two-dimensional digital filter 

dedicated to enhancing the profile is used to calculate 
the estimated value (hereinafter designated by a) of 
the noteworthy pixel A, the estimating section 56 
applies the multilevel values of the pixels covered 

20 by the matrix 100b of FIG. 4 (a, b, c, d, e, f, g, h, 
i in FIG. 4) to the two-dimensional digital filter 
dedicated to enhancing the profile (e.g., Laplacian 
filter or Prewitt filter), shown in FIGS. 5A and 5B, 
and FIGS. 6A and 6B. 

25 For example, when the Laplacian filter of FIG. 

5A is used, the estimated value a is given by 
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a=(axO)+{b*(-l) }+(cxO)+{dM-D }+ (ex4) + 
{fx(-l) }+(gxO)+{hx(-l)}+(ixO) . 



Also when using the Laplacian filter of FIG. 5B, 
5 the estimated value a can be obtained likewise. 

In the illustrated embodiment, two types of 
Laplacian filters of FIGS. 5A and 5B, are used. The 
present invention should by no means be limited to this 
illustrated examples, and alternative Laplacian 
10 filters having other patterns may be used. 

Specifically, as long as it satisfies all the following 
conditions (i) through (iv), any two-dimensional 
digital filter dedicated to enhancing the profile can 
give the same results as those with a Laplacian filter. 
15 Here, alphabetical characters a through i designate 
pixel values each corresponding to a respective pixel 
value in the matrix 100b of FIG. 4. 

(i) | e | = I a + b + c + d + f + g + h + i | 

(ii) The sign of e is opposite to that of Jb, d, 

20 £, h. 

(iii) The sign of e is the same as that of a, 
c, g, 1; ora=c=g=i=0. 

(iv) In the mat rix 10 Ob (FIG. 4), the values each 
obtained by multiplying the distance between the 

25 centering pixel (pixel value e) and each of the ad j acent 
pixels (pixel values a, Jb, c, d, £ , g, h, 1) by |a|, 
|b|, | c | , I d | , I f I , Igl, |h|, | i | , respectively, become 
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close to one another in largeness of digit (except that 
a, c r g , 1 = 0 ) . 

Therefore, by employing the Laplacian filter of 
FIGS. 5A and 5B or any of two-dimensional digital 
5 filters dedicated to enhancing the profile that satisfy 
all the foregoing conditions (i) through (iv) , it is 
possible to estimate the multilevel pixel value (the 
estimated value or) of the noteworthy pixel A. 

The manner in which the estimated value or is 
10 calculated using the Prewitt filters of FIGS. 6A and 
6B, will now be described. In an example that uses 
the Prewitt filter of FIG. 6A, the estimated value a 
is given by 



a=[[(axl)+(bxO)+{cx(-l)}+(dxl)+(exO)+ 
{fx(-l)}+(gxl)+(hxO)+{ix(-l)}] 2 +[(axl)+ 

(bxl)+( C xl)+(dxO)+(exO)+(fxO)+{gx(-l)}+ 
{hx (-1) } + {ix (-1) }] 2 ] 172 . 



20 Also in another example using the Prewitt filter 

of FIG. 6B, the estimated value a is obtained in the 
same way. 

The estimating section 56 outputs the operation 
result, or the estimated value or, which is obtained 
25 by applying the two-dimensional digital filter 

dedicated to enhancing the profile to the pixels in 
the predetermined area, as the pixel value of the 
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noteworthy pixel A. 

(3) Third Estimating Method: 

The estimating section 56 estimates the value 
of the noteworthy pixel A based on the values of 
5 non-binari zed pixels to which an error is diffused by 
the error diffusion section 52. 

FIGS. 7A and 7B and FIGS. 8A and 8B are diagrams 
for describing the manner in which the value of the 
noteworthy pixel A is estimated based on the values 
10 of non-binarized pixels to which the error is diffused. 

FIG. 7A illustrates the non-binarized pixels 
neighboring the noteworthy pixel A; FIG. 7B, a 
two-dimensional digital filter to be applied to the 
pixels of FIG. 7A. 
15 Assuming that the non-binarized pixels 

neighboring the noteworthy pixel A respectively have 
the multilevel values of a, b, c, d as shown in FIG. 
7A, the estimating section 56 applies the 
two-dimensional digital filter having weighting 
20 factors as of FIG. 7B to the pixel values a, b, c, d 
to estimate the noteworthy pixel A. 

Specifically, in the example of FIGS. 7A and 7B, 
the multilevel value a of the noteworthy pixel A is 
calculated with the following equation (i) as the 
25 estimated value a of the noteworthy pixel A. 

a=(a + b/ 2 + c + d/ 2)/3 . . . (i) 
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The two-dimensional digital filter of FIG. 7B 
is preferred when performing an error diffusion 
technique {FIG. 10A) well known as "Floyd and 
5 Steinberg's Method". 

FIG. 8A illustrates the non-binar i zed pixels 
adjacent to the noteworthy pixel A; FIG. 8B, a 
two-dimensional digital filter to be applied to the 
pixels of FIG. 8A. 
10 Assuming that the non-binar i zed pixels adjacent 

to the noteworthy pixel A respectively have the 
multilevel values of a, b, c, d, &, f, g, h, ±, j r k, 
1 as shown in FIG. 8A, the estimating section 56 applies 
the two-dimensional digital filter having weighting 
15 factors as of FIG. 8B to the pixel values a through 
1 to estimate the noteworthy pixel A. 

In FIG. 8B, the weighting factors are designated 
by reference characters X, Y, Z, U: for example, X = \ , 
Y = l/2, Z = l/V2, C/ = l/V3. These weighting factors X, 
20 Y, Z, U should by no means be limited to these values 
and may be changed to other values. 

In the example of FIGS. 8A and 8B, the multilevel 
value a of the noteworthy pixel A is calculated with 
the following equation (ii) to serve as the estimated 
25 value a of the noteworthy pixel A. 
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a = {X-(a + e)+Y-(b + j) + Z-(h + l) / 2 + U -(c + g + i + k)} 

^{(l8 + 9x2 1/2 +8x3 1/2 ) / 6} '" Ul) 

The two-dimensional digital filter of FIG. 8B 
is preferred when performing error diffusion 
5 techniques ( described later with reference made to FIGS . 
11A and 11B) well known as "Javis, Judice and Ninkes' 
Method" and "Stucki's Method" (described later with 
reference made to FIGS. 11A and 11B). 

In the present embodiment , the estimating section 

10 56 uses one of the above three types of methods (1) 
through (3) to calculate the estimated value a of the 
noteworthy pixel A. 

Thebinarizing section 51 converts the multilevel 
value of a noteworthy pixel currently under scanning 

15 into a binary value based on the estimated value a of 
the noteworthy pixel A that is calculated by the 
estimating section 56. Precisely, the estimated 
value a of the noteworthy pixel A thus calculated by 
the estimating section 56 is compared with a 

20 predeterminedvalue (e.g., 128) ; ifthe estimated value 
a of the noteworthy pixel A is equal to or larger than 
the predetermined value (<x^ 128), 255 (ON) is recorded 
as the value of the noteworthy pixelAina predetermined 
area of the RAM 16 or the magnetic disc 26. Hereinafter, 

25 thepixelvalueofthe noteworthy pixel A obtained after 
binarizing is designated by at' (cc' = 0 or a' = 255) . 
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To the contrary, if the pixel value a of the 
noteworthy pixel A is smaller than the predetermined 
value (a<128) , 0 (OFF) is recorded as the pixel value 
a' of the noteworthy pixel in a designated area of the 
5 RAM 16 or the magnetic disc 26. 

The pixel-on-prof ile detection section 54 
discriminates whether or not the noteworthy pixel is 
a pixel that constitutes part of a profile of the 
multilevel input image (for convenience to describe, 
10 this discrimination will hereinafter occasionally be 
called "profile detection") . The pixel-on-prof ile 
detection section 5 4 calculates a profile value E with 
respect to the noteworthy pixel A based on both the 
multilevel value of the noteworthy pixel and those of 
15 the adjacent pixels, using a Laplacian filter or a 
Prewitt filter or by making addition and subtraction 
individually on the multilevel values of the noteworthy 
pixel A and the adjacent pixels. 

The thus obtained profile valueEis then compared 
20 with a predetermined value, or a threshold T , to 

discriminate whether or not the noteworthy pixel is 
the pixel that constitutes part of the profile of the 
multilevel input image. 

The threshold T, which can be set by an operator, 
25 plays as a parameter to control the degree of 
enhancement of the profile. 

When calculating the profile value Ewith respect 
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to the noteworthy pixel A, the pixel-on-prof ile 
detection section 54 scans the multilevel image data 
expanded on a memory ( RAM 1 6 ) while arranging thematrix 
100b (FIG. 4) so as to center the noteworthy pixel A 
5 (pixel value e of FIG. 4) . For every noteworthy pixel 
A centering the matrix 100b, the profile value E is 
calculated using the digital filter dedicated to 
enhancing the profile. The profile value E is 
calculated in the same way as the estimated value a. 

10 When the two-dimensional digital filter 

dedicated to enhancing the profile is used to calculate 
t he profile value Ewith respect to the noteworthy pixel 
A, both the multilevel value (pixel value e) of the 
noteworthy pixel A and those of the adjacent pixels 

15 in the matrix 100b of FIG. 4 (a, b, c, d, f, g, h, i 
in FIG. 4) are applied to the two-dimensional digital 
filter dedicated to enhancing the profile (e.g., 
Laplacian filter or Prewitt filter) , which filter is 
shown in FIGS. 5A and 5B, and FIGS. 6A and 6B. 

20 For example, when the Laplacian filter of FIG. 

5A is used, the profile value E is given by 

E=(a*0)+{b*(-l) }+(cxO)+{dx(-l) }+(e*4)+ 
{fx(-l) }+(gx0)+{hx(-l) }+(ix0) . 

25 

Also when using the Laplacian filter of FIG. 5B, 
the profile value E can be obtained likewise. 
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When the pixel-on-prof ile detection section 54 
calculates the profile value E, two types of Laplacian 
filters of FIGS. 5A and 5B are used. The present 
invention should by no means be limited to this 
5 illustrated examples, and alternative Laplacian 

filters having other patterns may be used. Precisely, 
as long as it satisfies all the following conditions 
(i) through (iv), any two-dimensional digital filter 
dedicated to enhancing the profile can give the same 
10 results as those with a Laplacian filter. Here, 

alphabetical characters a through i designate pixel 
values each corresponding to a respective pixel value 
in the matrix 100b of FIG. 4. 

15 (i) | e | = I a + b + c + d + f + g + h + i I 

(ii) The sign of e is opposite to that of i>, d, 

£ , h. 

(iii) The sign of e is the same as that of a, 
c, g, 1; or a=c=g=i=0. 

20 (iv) In the matrix 100 (FIG. 4), the value s each 

obtained by multiplying the distance between the 
noteworthy pixel (the pixel value e) and each of the 
adjacent pixels (pixel values a, h, c, d, £, g, h, 1) 
by |a|, |b|, |c|, |d|, |f|, |g|, |h|, |i|, respectively, 

25 become close to one another in largeness of digit 
(except that a, c, g, 1=0). 
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Therefore, by employing the Laplacian filter of 
FIGS. 5A and 5B, or any of two-dimensional digital 
filters dedicated to enhancing the profile that satisfy 
all the foregoing conditions (i) through (iv) , it is 
5 possible to discriminate whether or not the noteworthy 
pixel is a pixel on the profile of the image. 

The manner in which the profile value E is 
calculated using the Prewitt filters of FIGS. 6A and 
6B will now be described. In an example using the 
10 Prewitt filter of FIG. 6A, the profile value E is given 
by 

E=[[(axl)+(bxO)+{cx(-l)}+(dxl)+(exO)+ 
{fx(-l) }+(gxl) + (hxO)+{ix(-l) }] 2 +[(axl) + 
15 (bxl)+( C xl)+(dxO)+(exO)+(fxO)+{gx(-l)}+ 
{hx (-1) } + {ix (-1) } ] 2 ] 1/2 . 

Also in another example using the Prewitt filter 
of FIG. 6B, the profile value E is obtained in the same 
20 way. 

The greater the variation among the pixel values 
of the noteworthy pixel and the adjacent pixels, the 
greater the profile value E will be obtained. 

Then the pixel-on-prof ile detection section 54 
25 compares the thus obtained profile value E with the 
predetermined threshold T . If the profile value E is 
greater than the threshold T , it is judged that the 
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noteworthy pixel (the pixel value e) is a pixel on the 
profile of the image; this is, the profile has been 
detected . 

The direction-of -prof ile detection section 55 
detects the direction in which the profile of the 
multilevel input image extends with respect to the 
noteworthy pixel A. The manner in which the 
direction-of -prof ile detection section 55 detects the 
direction of the profile will be described with 
reference to FIG. 9, which illustrates various profile 
directions . 

The direction-of -prof ile detection section 55 
makes the following calculations using the pixel values 
a, b, c, d, e, f, g, h, 1 in the matrix 100b (see FIG. 
4), which centers the noteworthy pixel. 

Xi = b + e + h 
X 2 = c + e + g 
X 3 = a + e + i 
X 4 = d + e + f 

The thus obtained values Xi, X 2 , X 3 , X 4 are compared 
with one another, and j (j is a natural number 
representing any of 1, 2, 3, 4) satisfying the following 
equation : 

Xj = m ax(Xi, X 2 , X 3 , X 4 ) 
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is calculated. If two or more of the values Xi, X 2 , 
X 3 , X 4 are equal to one another, the largest number 
j is adopted. 

5 Based on the thus obtained j, the direct ion-of - 

profile detection section 55 determines the profile 
direction Dj (Di, D 2 , D 3 , D 4 ) passing over the noteworthy 
pixel A (pixel value e) as shown in FIG. 9. 
Specifically in the matrix 100b, the direction passing 

10 along a string of the pixels b, e, h is designated as 
the profile direction Di; the direction passing along 
a string of the pixels c, e, q, the profile direction 
D 2 ; the direction passing along a string of the pixels 
a, e, i, the profile direction D 3 ; the direction passing 

15 along a string of the pixels d, e, £, the profile 
direction D 4 . From these profile directions Di, D 2 , 
D 3/ D4, a particular profile direction Dj is selected. 

In the halftoning apparatus of the present 
embodiment, the two-dimensional digital filter 

20 dedicated to enhancing the profile is used to calculate 
the profile value E. The present invention should by 
no means be limited to this illustrated example, and 
various changes or modifications may be suggested 
without departing from the gist of the invention. 

25 The technique will now be described in which the 

profile value E is calculated without the 
two-dimensional digital filter dedicated to enhancing 
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the profile, simultaneously with detecting the profile 
direction. The alphabetical characters a through i 
designate the respective pixel values; like reference 
numbers designate similar parts or elements throughout 
5 several views , so their detailed description is omitted 
here . 

(i) In the 3x3 matrix 100b of FIG. 4, the 
following calculations ( addition and sub t r act ion ) are 
made : 

10 
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20 (ii) Xi through X 8 are examined one by one to 

find the maximum of them. 

(iii) Based on the result of (ii), the profile 
direction is determined as follows: 

a) If Xi or X2 is the maximum, the location 
25 and the direction of the profile are Di (FIG. 

9) - 

b) If X3 or X 4 is the maximum, the location 
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and the direction of the profileareD 2 (FIG. 
9) . 

c) If X 5 or X 6 is the maximum, the location 
and the direction of the profileare D 3 (FIG. 
9) . 

d) If X 7 or X 8 is the maximum, the location 
and the direction of the profile are D 4 (FIG. 
9) . 

(iv) The location and the direction of the 
profile at the time of detection of the maximum are 
stored in the RAM 16, the magnetic disc 26, or the like. 

(v) The maximum of Xi through X 8 obtained in (ii) 
above is multiplied by 256/768 ( = 1/3). Then the re suit 
is rounded to the nearest whole number and stored in 
the RAM 16, the magnetic disc 26, or the like. 

In this manner, the profile value E can be 
calculated by making addition and subtraction 
individually on the values of the noteworthy pixel and 
the adjacent pixels without the two-dimensional 
digital filter dedicated to enhancing the profile. 
Alternatively, the profile value E may be calculated 
without making any multiplication or division so that 
the load on the CPU 10 can be decreased, speeding up 
calculation of the profile value E and, at the same 
time, detecting the profile direction. 

The value "256/768" by which the maximum of Xi 
through X 8 is multiplied in (v) above will now be 
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described . 

In order to obtain the value X ± (e.g., Xi) , the 
value obtained by adding the values of three pixels 
is subtracted from that obtained by adding the values 
5 of the remaining three pixel values 

(a + d + g- b- e-h). As the pixel values a through 
i are in the range of from 0 to 255 (8 bits) , the maximum 
and the minimum of X± (a + d+g-b-e-h) are 765 
(= 255 x 3) and 0, respectively. 
10 To control the value X± using the 8-bit threshold 

T, the number by which X ± is to be multiplied takes 
256 /(256 x 3 = 768 ) as a relatively small normalizing 
factor . 

The error diffusing section 52 diffuses an error 
15 Z, which has occurred in the binary value with respect 
to the noteworthy pixel as the binarizing section 51 
performs the binarizing, to unscanned pixels adjacent 
to the noteworthy pixel, 
where , 

20 Z = 255-e (if e^l28) 

Z=e (if e<128) 

(e: the pixel value of the noteworthy pixel). 

More specifically, the error diffusing section 
25 52 diffuses the error Zto the unscanned pixels adjacent 
to the noteworthy pixel A by one diffusion technique. 
The occurred error Z with respect to the noteworthy 
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pixel A is proportionally distributed to the plural 
unscanned pixels adjacent to the noteworthy pixel e 
in accordance with a weighting pattern (weighting 
factors) as of FIGS. 10 and 11. 
5 FIGS. 10A through 10D, show the weighting 

patterns of the error diffusion technique: in what 
ratio the error Z with respect to the noteworthy pixel 
A, which is designated by • in FIGS. 10A through 10D, 
is distributed to the unscanned pixels adjacent to the 
10 noteworthy pixel A. 

The error diffusing based on the weighting 
pattern of FIG . 1 OA will now be described . For example, 
if the pixel value e of the noteworthy pixel A is e = 100 
(e<128), the error Z=100; consequently 100x7/16 is 
15 added to an unscanned pixel (pixel value f in FIG. 4) 
neighboring the noteworthy pixel A in the primary 
scanning direction. 

Likewise, 100x5/16 is added to another unscanned 
pixel (pixel value h in FIG. 4) neighboring the 
20 noteworthy pixel A in the secondary scanning direction 
(perpendicular to the primary scanning direction) ; 
100x1/16, still another unscanned pixel (pixel value 
i in FIG. 4) neighboring the noteworthy pixel A in the 
primary scanning direction; 100x3/16, a further 
25 unscanned pixel (the pixel value g in FIG. 4) 

neighboring the noteworthy pixel A in the primary 
scanning direction. As a result, the occurred error 
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is proportionally distributed (diffused) to the 
unscanned pixels adjacent to the noteworthy pixel A. 

In addition, the error diffusing section 52 has 
also other weighting patterns as of FIGS. 10B through 
5 10D, so that the error Z can be proportionally 

distributed to the unnscanned pixels adjacent to the 
noteworthy pixel A based on these plural (4 in the 
present embodiment) types of weighting patterns , thus 
realizing the error diffusing. 
10 Also when using the weighting patterns of FIGS. 

10B through 10D, the error diffusing section 52 
performs the error diffusing in the same manner as when 
using the weighting pattern of FIG. 10A. 

Further, when the pixel-on-prof ile detection 
15 section 54 detects that the noteworthy pixel A is a 
pixel constituting part of a profile of the multilevel 
input image, the error diffusing section 52 may perform 
an exceptional process as described below. 

As this exceptional process, values according 
20 to the occurred error Z are added to the values of the 
unscanned pixels along the profile in the detected 
direction, which is detected by the 
direction-of -prof ile detection section 54. 
Specifically, the error diffusing section 52 decides 
25 the profile direction Dj , which has been detected by 
the direction-of -prof ile detection section 54, and 
then performs the following processes in accordance 
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with the decision. 

(i) If the profile direction is Di, Z/2 is added 
to the pixel value h. 

(ii) If the profile direction is D 2 , Z/2 is added 
5 to the pixel value g. 

(iii) If the profile direction is D 3 , Z/2isadded 
to the pixel value 1 . 

(iv) If the prof ile direction is D 4 , Z/2 is added 
to the pixel value f. 

10 

In the present embodiment, the value obtained 
by dividing the error Z by 2 is added to the pixel value 
of the unscanned pixel along the profile direction. 
The present invention should by no means be limited 

15 to this, and various changes or modifications may be 
suggested without departing from the gist of the 
invention. As an alternative example, the value 
obtained by dividing the error Z by a number other than 
2 may be added to the pixel value of the unscanned pixel 

20 along the profile in the detected direction in 

accordance with printing conditions such as the quality 
of ink (toner) or a printing medium to be used. 

Meanwhile, binarizing is realized without error 
diffusing by the error diffusion section 52 in the 

25 half toning apparatus of the present embodiment . Such 
binarizing will be hereinafter called "simple 
threshold (binarizing) method". 
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The error diffusion technique changing section 
53 changes the technique of the error diffusing, which 
is currently used by the error diffusing section 52, 
to another in accordance with a predetermined way as 
5 the scanning of thepixels oft he multilevel input image 
progresses. For example, the plural weighting 
patterns of FIGS. 10A through 10D are selectively 
changed from one to another in the sequence of A, B, 
C, D, A, B. . . . Namely, the error diffusion technique 
10 changing section 53 makes such change by changing the 
weighting pattern currently used in the error diffusing 
section 52 to another that is selected in the 
predetermined sequence. Alternatively the weighting 
patterns may be changed at random. 
15 Furthermore, the error diffusion technique 

changing section 53 changes the error diffusion 
technique to another when the pixel-on-prof ile 
detection section 54 judges that the noteworthy pixel 
is apixel constituting part of theprofile of the image 
20 (around the profile) . Otherwise the error diffusion 
technique changing section 53 changes the error 
diffusion technique for every pixel. 

The way of change of error diffusion technique 
by the error diffusion technique changing section 53 
25 is exemplified by (1) changing the weighting pattern, 
and (2) changing the error diffusion technique it self ; 
either changing way may be selected according to the 
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need. In the present embodiment, the changing of the 
error diffusion technique by the error diffusion 
technique changing section 53 can be omitted to perform 
the error diffusing with one particular technique. 
5 The operation of the individual element or part 

of the halftoning apparatus of the present embodiment 
during the change of the error diffusion technique 
(weighting pattern) around the profile will now be 
described with reference to the flowchart of FIG. 12, 

10 steps A10 through A110. 

First of all, the halftoning apparatus binarizes 
the multilevel image data expanded on the RAM 16 while 
scanning the image data using the matrix 100b, which 
centers the noteworthy pixel A, as shown in FIG. 4. 

15 At that time, in the matrix 100b, the pixel-on-prof ile 
detection section 54 first calculates the profile value 
E with respect to the noteworthy pixel A based on both 
the multilevel value of the noteworthy pixel A and those 
of the adjacent pixels, by using the two-dimensional 

20 digital filter dedicated to enhancing the prof ile , such 
as the Laplacian filters of FIGS. 5A and 5B, or the 
Prewitt filters of FIGS . 6Aand6B, or by ma king addi t i on 
and subtraction individually on the multilevel values 
of the noteworthy pixel A and the adjacent pixels for 

25 direct calculation (step A10). 

Then, the binarizing section 51 compares the 
estimated value or of the noteworthy pixel A, which is 



38 



calculated by the estimating section 56, with the 
predetermined value (128 in the present embodiment) 
(step A20) . If the estimated value a is equal to or 
larger than 128 (a ^128) (the YES route of step A20) , 
5 255 is recorded as the pixel value <x r of the noteworthy 
pixel A in a predetermined area of the RAM 16 or the 
magnetic disc 26 (step A30). 

And, 255 minus a is calculated. The result of 
this calculation is recorded as the error Z in the 

10 predetermined area of the RAM 16 or the magnetic disc 
26 (step A40) . 

Otherwise if the estimated value a is smaller 
than 128 (e<128) (the NO route of step A20), the 
binarizing section 51 records 0 as the pixel value a' 

15 of the noteworthy pixel A in the predetermined area 
of the RAM 16 or the magnetic disc 26 (step A50) ; and 
the binarizing section 51 records the estimated value 
a of the noteworthy pixel A as the error Z in the 
predetermined area of the RAM 16 or the magnetic disc 

20 26 (step A60) . 

Next the pixel-on-prof ile detection section 54 
compares the profile value E with a preset threshold 
T (step A70) . If the profile value E is smaller than 
the threshold T (the NO route of the step A70), the 

25 error diffusing section 52 proportionally distributes 
the error Z to the plural un scanned pixels in accordance 
with the weighting pattern of the current error 
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diffusion technique (step A80). 

On the contrary, if the profile value E is equal 
to or larger than T (see the YES route of step A70), 
the exception process is carried out (step A100); 
5 namely, the profile direction Dj is detected by the 
direct ion-of -prof ile detection section 55, and then 
the values (Z/2) corresponding to the occurred error 
are added to the pixel values of the unscanned pixels 
along the profile in the detected direction. 

10 After that, the error diffusion technique 

changing section 53 selectively changes the weighting 
pattern from one to another among the plural weighting 
patterns of FIGS. 10A through 10D, (for example, in 
the sequence of A, B, C, D, A, B . . . in FIG. 10) (step 

15 A110) . 

After completion of step A8 0 or A110, for the 
multilevel image data expanded on the RAM 16, the 
halftoning apparatus shifts the matrix 100b so as to 
center the next noteworthy pixel (the pixel next to 
20 the current noteworthy pixel in the scanning direction) 
(step A90) , and then returns to step A10 to repeat the 
same processes for a newly given noteworthy pixel. 

In the present embodiment, as mentioned above, 
the change of the error diffusion technique by the error 
25 diffusion technique changing section 53 is made by 
selecting one from the plural weighting patterns of 
FIGS. 10A through 10D, in the sequence of A, B, C, D, 
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A, B.... Alternatively selection from the plural 
weighting patterns of FIGS. 10A through 10D, may be 
at random. 

In the halftoning apparatus of the present 
5 embodiment, the error diffusion technique changing 
section 53 changes the error diffusion technique only 
when the pixel-on-prof ile detection section 54 judges 
that the noteworthy pixel A is a pixel on the profile 
of the input multilevel image. Alternatively the 

10 change of the error diffusion technique may be made 
for every pixel. In this case also, a desired error 
diffusion technique (weighting pattern) may be 
selected from the plural candidates at random. 

The operation of the individual element or part 

15 of the halftoning apparatus during the change of the 
error diffusion technique (weighting pattern) for 
every pixel will now be described with reference to 
the flowchart of FIG. 13, steps BIO through B110. 

First of all, the halftoning apparatus binarizes 

20 the multilevel image data expanded on the RAM 16 while 
scanning the image data using the matrix 100b of FIG. 
4, which centers the noteworthy pixel. At that time, 
in the matrix 100b, the pixel-on-prof ile detection 
section 54 first calculates the profile value E with 

25 respect to the noteworthy pixel based on both the 
multilevel value of the noteworthy pixel and those of 
the adjacent pixels, by using the two-dimensional 
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digital filter dedicated to enhancing the profile, such 
as the Laplacian filter (FIGS. 5Aand5B) or the Prewitt 
filter (FIGS. 6A and 6B), or by making addition and 
subtraction individually on the multilevel values of 
5 the noteworthy pixel A and the adjacent pixels for 
direct calculation (step BIO). 

Then the binarizing section 51 compares the 
estimated value a of the noteworthy pixel A, which is 
calculated by the estimating section 56, with the 

10 predetermined value (128 in the present embodiment) 
(step B20) . If the estimated value a is equal to or 
larger than 128 (a^l28) (the YES route of step B20), 
255 is recorded as the pixel value a' of the noteworthy 
pixel A in a predetermined area of the RAM 16 or the 

15 magnetic disc 26 (step B30) . 

And, 255 minus a is calculated, and the result 
of this calculation is recorded as the error Z in the 
predetermined area of the RAM 16 or the magnetic disc 
26 (step B40) . 

20 Otherwise if the pixel value or is smaller than 

128 (a<128) (the NO route of step B20) , the binarizing 
section 51 records 0 as the pixel value a f of the 
noteworthy pixel A in the predetermined area of the 
RAM 16 or the magnetic disc 26 (step B50); and the 

25 binarizing section 51 records the estimated value a 
of the noteworthy pixel as the error Z in the 
predetermined area of the RAM 16 or the magnetic disc 
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26 (step B60) . 

Next the pixel-on-prof ile detection section 54 
compares the profile value E with a preset threshold 
T (step B70) . If the profile value E is smaller than 
5 the threshold T (the NO route of the step B70), the 
error diffusing section 52 proportionally distributes 
the error Z to the plural unscanned pixels adjacent 
to the noteworthy pixel A in accordance with the 
weighting pattern of the current error diffusion 
10 technique (step B80) . 

After that, the error diffusion technique 
changing section 53 selectively changes the plural 
weighting patterns of FIGS. 10A through 10D, one after 
another (e.g., in the sequence of A, B, C, D, A, B...) 
15 (step B90) . 

On the contrary, if the profile value E is equal 
to or larger than T (see the YES route of step B70), 
the exception process is carried out (step B110); 
namely, the profile direction Dj is detected by the 
20 direction-of -prof ile detection section 55 and then the 
values (Z/2) corresponding to the occurred error are 
added to the pixel values of the unscanned pixels along 
the profile in the detected direction. 

After completion of step B90 or B110, for the 
25 multilevel image data expanded on the RAM 16, the 
halftoning apparatus shifts the matrix 100b so as to 
center the next noteworthy pixel A (step B100), and 
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then returns to step BIO to repeat the same processes 
for a newly given noteworthy pixel A. 

As also mentioned above, the error diffusion 
technique ( t he we ight ing pa 1 1 er n ) is changed for every 
5 pixel, thus minimizing moire and artifacts in the 
binary image, which did not appear in the original 
(multilevel) image. 

How to change the current error diffusion 
technique to another will now be described with 

10 reference to FIGS . 11A and 11B. FIG. 1 1A illustrates 
the weighting pattern of the error diffusion technique 
by Javis, Judice and Ninke, while FIG. 11B illustrates 
the error diffusion technique by Stucki. 

In these error diffusion techniques of FIGS. 11A 

15 and 11B, error Z proportionally distributes to the 
noteworthy pixel and the unscanned pixels adjacent to 
the noteworthy pixels. In these techniques the area 
ofpixelstowhichtheerrorZis distributed is enlarged 
as compared to the error diffusion technique of FIG. 

20 10 . 

The weighting pattern of FIG. 11A is known as 
the error diffusion technique of Javis, Judice and 
Ninkes, while that of FIG. 11B is known as the error 
diffusion technique of Stucki. 
25 In these techniques , t he e r r or di f f us ing s ect ion 

52 proportionally distributes (diffuses) the error Z 
using either the weighting pattern of FIG. 11A or that 
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of FIG . 11B. And the error diffusion technique 
changing section 53 changes the weighting patterns one 
another. The remaining processes are substantially 
identical with those illustrated in the flowchart of 
5 FIG. 12, steps A10 through A110. 

With this arrangement, scanning the pixels 
contained in the multilevel input image, the halftoning 
apparatus of the present embodiment converts 
(binarizes) the multilevel pixel values of the image 

10 data into the binary values for every noteworthy pixel 
A being currently scanned. 

Here, the profile detection process and the 
binarizing process carried out in the halftoning 
apparatus of the present embodiment will now be 

15 described with reference made to the flowchart (step 
CIO to C120) of FIG. 14. 

First of all, the operator inputs on the operation 
panel a threshold T as a preset parameter controlling 
the degree of profile enhancement ( step CIO ) , whereupon 

20 the halftoning apparatus starts raster scanning (step 
C20) . 

At that time, the binarizing may start from any 
of four pixels respectively located at four corners 
of the input image, and may progress in a desired 
25 direction. 

Then discrimination is made whether or not the 
scanning has been completed (step C30) ; if the result 
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of this discrimination is positive (the YES route of 
step C30) the operation of the halftoning apparatus 
is terminated. 

Otherwise if the scanning is to be continued (the 
5 NO route of step C30) , the halftoning apparatus selects 
one particular way from the plural profile detection 
ways; i.e. (1) using the two-dimensional digital filter 
dedicated to enhancing the profile, and (2) directly 
calculating from the pixel values (step C40). 

10 Alternatively the profile detect ion way may be selected 
by the operator. 

Here, if a Laplacian filter (the two-dimensional 
digital filter dedicated to enhancing the profile) is 
selected to use, the profile value E is calculated using 

15 the Laplacian filter as mention above, and the obtained 
profile value E is compared with the threshold T for 
judgment that the noteworthy pixel A being currently 
scanned is a pixel constituting part of the profile 
of the image (step C50) . 

20 Likewise, if a Prewitt filter (the 

two-dimensional digital filter dedicated to enhancing 
the profile) is selected to use, the profile value E 
is calculated using the Prewitt filter as mentioned 
above, and the obtained profile value E is compared 

25 with the threshold T for judgment that the noteworthy 
pixel A being currently scanned is a pixel constituting 
part of the profile of the image (step C60). 
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When the direct calculation from the pixel values 
is selected, the profile value Eis calculated by ma king 
addition and subtraction individually on the 
multilevel values of the noteworthy pixel A and the 
5 adjacent pixels, and then the obtained profile value 
E is compared with the preset threshold T for judgment 
that the noteworthy pixel A being currently scanned 
is a pixel constituting part of the profile of the image 
(step C70) . 

10 Next, the half toning apparatus selects the method 

of halftoning (stepC80). Specifically, aparticular 
one is selected from the following three halftoning 
methods: with plural weighting patterns of the error 
diffusion technique provided in advance, (1) changing 

15 the weighting pattern only upon detection that the 
noteworthy pixel A being currently scanned is a pixel 
on the profile of the image (only upon encountering 
of the profile); (2) changing the weighting patterns 
one to another for every pixel; and (3) changing the 

20 weighting pattern to that of another error diffusion 
technique as shown in FIGS. 11A and 11B. 

If the method (1) above is selected (step C90), 
one identical weighting pattern is constantly used in 
halftoning the pixels lying between one profile and 

25 the other prof ile . At that time , the weight ing pattern 
is preferably changed either in a predetermined 
sequence or at random. 
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Otherwise if the method (2) above is selected 
(step C100) , the weighting pattern is changed for every 
pixel; at that time, the weighting pattern is 
preferably changed either in a predetermined sequence 
5 or at random. 

Still otherwise if the method (3) is selected, 
the halftoning is performed using such error diffusion 
methods as shown in FIGS. 11A and 11B (step C110). 
After that, for the multilevel image data 
10 expanded on the RAM 16, the halftoning apparatus shifts 
the matrix 100 so as to center the next noteworthy pixel 
(step C120), and then returns to step C30. 

Subsequently the processes (steps CIO to C120) 
are repeated until the raster scanning is completed. 
15 Next, the entire operation, including pixel value 

estimation, of the halftoning apparatus of the present 
embodiment will now be described with reference made 
to the flow chart (steps D10 to D120) of FIG. 15. 

First of all, an operator inputs a weighting 
20 factor and an area to be used for estimating the 

multilevel value of the noteworthy pixel A through a 
keyboard 20 or a mouse 22 (step D10), whereupon the 
halftoning apparatus starts raster scanning (step 
D20) . 

25 At that time, binarizing process may start from 

any of four pixels respectively located at four corners 
of the input image, and may progress in a desired 
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direction . 

Then discrimination is made whether or not the 
scanning has been completed (step D30) ; if the result 
of this discrimination is positive (the YES route of 
5 step D30) the operation of the halftoning apparatus 
is terminated. 

Otherwise if the scanning is to be continued (the 
NO route of step D30) , the halftoning apparatus selects 
one particular way from the above-mentioned three ways 

10 of calculating the estimated value a (step D40). 

Alternatively the calculation way may be selected by 
the operator. 

Here, upon selection of the use of a 
two-dimensional digital filter other than the digital 

15 filter dedicated to enhancing the profile, the 
estimated value a is calculated with the 
two-dimensional digital filter in the same manner as 
described above under the caption of "First Estimating 
Method" (step D50) . 

20 Otherwise if the two-dimensional digital filter 

dedicated to enhancing the profile is selected to use, 
the estimated va lue a is calculated using the Lap lacian 
filter or the Prewitt filter in the same manner as 
described above under the caption of "Second Estimating 

25 Method" (step D60) . 

Still otherwise if selected the calculation based 
on the values of non-binarized pixels to which the error 
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is diffused, the estimated value a is calculated, 
applying the two-dimensional digital filter to those 
non-binarized pixels, in the same manner as described 
above under the caption of "Third Estimating Method" 
5 (step D70) . 

Next, the halftoning apparatus selects the method 
of halftoning (step D80) . Specifically, a particular 
one is selected from the following three types of 
halftoning methods: (1) the simple threshold method; 
10 (2) the method in which a common error diffusion 

technique is used; (3) the method in which the profile 
detection (shown in FIG. 14) is per formed before error 
diffusing . 

Upon selection of the above method (1) (step D90) , 
15 the binarizing section 51 binarizes the estimated 
multilevel value a, which is calculated by the 
estimating section 56 as the pixel value of the 
noteworthy pixel A. At that time, the error {the 
difference between the po s t -binar i z ing pixel value ct r 
20 (= 0 or 255) and the estimated value a} accompanying 
the binarizing is not diffused (simple threshold 
method) . 

Otherwise if the above method (2) is selected 
(step D100), the binarizing section 51 binarizes in 
25 the simple threshold method the estimated value or, 
which is calculated by estimating section 56 as the 
pixel value of the noteworthy pixel A, whereupon the 
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error diffusion section 52 diffuses the error {the 
difference between the pixel value ct' (= 0 or 255) after 
binarized and the estimated value or) accompanying the 
binarizing by using a well-known error diffusion 
5 technique. 

Still otherwise if the above method (3) is 
selected (step DUO), the binarizing section 51 
binarizes in the simple threshold method the estimated 
value a, which is calculated by estimating section 56 
10 as the pixel value of the noteworthy pixel A, whereupon 
steps C40 to C110 of the flowchart of FIG. 14 are 
performed. In this case, "weighting factor and area 
used to estimate the noteworthy pixel A" has been input 
in advance. 

15 After that, for the multilevel image data 

expanded on the RAM 16, the halftoning apparatus shifts 
the matrix 100b so as to center the next noteworthy 
pixel A (step D120), and then returns to step D30. 
Subseguently the processes (step D10 through 

20 dl20) are repeated until the raster scanning is 
completed . 

As described above, when converting the 
multilevel value of the note worthy pixel A into a binary 
value, the value to be converted is not the noteworthy 
25 pixel A itself but the estimated pixel value ct, which 
is calculated based on other pixel values than that 
of the noteworthy pixel A. It is thus possible to 
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soften the sudden shift in the pixel values even in 
a patterned image, such as a checkered one, in which 
neighboring two pixels can greatly differ in pixel 
values . 

5 The softening of the sudden shift in the pixel 

values reduces bands (moire, etc.) that appear in a 
halftoned image along the scanning direction, which 
bands are due to the sudden shift in the pixel values 
in an original image. 

10 Further, even if some pixels have outstanding 

pixel values due to image defects, significant errors 
risen from the outstanding value are prevented from 
being directly diffused to unprocessed adjacent pixels 
It is thus possible to minimize moire and other 

15 artifacts that did not appear in the original 

(multilevel) image. Also, even if the original image 
contains latent, not visually apparent, moire, which 
is due to the periodic structure typical of an input 
device (scanner, digital camera, etc) , it is prevented 

20 that the moire is enhanced by error diffusing to 

resultantly become visually apparent in the halftoned 
image. Moreover, even in a very gentle tone-gradation 
image, a peculiar artifact that tends to appear in such 
an image can be successfully prevented. 

25 The use of a two-dimensional digital filter 

facilitates the calculation of the estimated value a. 
In particular, a two-dimensional digital filter 
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dedicated to enhancing the profile, such as the 
Lapracian filter or the Prewitt filter, provides a 
sharp, profile-enhanced appearance in the processed 
image, thereby producing a visually appealing image. 
5 Further, since the calculation of the estimated 

value a is performed by the estimating section 56 in 
parallel with the binarizing process by the binarizing 
section 51, it is possible to reduce the processing 
time . 

10 Furthermore, when the pixel-on-prof ile 

detection section 54 detects that the noteworthy pixel 
is contained in the profile of the input image, the 
error diffusion technique changing section 53 changes 
the error diffusion technique (weighting pattern) to 

15 another. It is therefore possible to localize the 
undesirable effects {moire and artifacts) , which have 
been presumably caused by changing the error diffusion 
technique, to dissolve into the profile, thus 
minimizing moire and artifacts in the binary image, 

20 which did not appear in the original (multilevel) image. 
Further, when the pixel-on-prof ile detection section 
54 detects that the noteworthy pixel is contained in 
the profile of the input image, the error diffusing 
section 52 performs the exception process of adding 

25 the values (Z/2), in accordance with the error occurred 
with the binarizing of the noteworthy pixel, to the 
values of the unscanned pixels along the profile in 
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the direction detected by the direct ion-or-prof ile 
detection section 55. It is accordingly possible to 
minimize moire and artifacts in the binary image, which 
did not appear in the original (multilevel) image. 
5 Still further, discrimination is made whether 

or not the noteworthy pixel is a pixel on the profile, 
by comparing the profile value E with respect to the 
noteworthy pixel, which is calculated based on the 
values of the noteworthy pixel and its adjacent pixels, 

10 to the threshold T (e.g., 128). Accordingly it is 
possible not only to facilitate detecting the profile 
but also to enable detecting the profile while scanning 
the multilevel image, realizing a high-speed profile 
detecting process. 

15 Further, with the halftoning apparatus of the 

present embodiment, the operator can change the degree 
of prof ile enhancement (sharpness) in the binary image 
by changing the threshold T, thereby obtaining printed 
matter devoid of a blur, which would have been caused 

20 in the conventional halftoning. 

For example, the closer the threshold T to 255, 
the more the profile will be enhanced; the closer the 
threshold T to 0, the more the profile will be blurred. 
Furthermore, the pixel-on-prof ile detection 

25 section 54 employs the two-dimensional filter 

dedicated to enhancing the profile, such as a Laplacian 
filter or a Prewitt filter to calculate the profile 
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value, thus facilitating detection of the profile. 

In addition, as the error diffusing section 52 
proportionally distributes the occurred error to the 
plural (4 in the present embodiment) unscanned pixels 
5 adj a cent to the noteworthy pixel basedonthe weighting 
patterns as shown in FIGS. 7A through 7D, the error 
diffusion technique changing section 53 selects one 
from the candidate weighting patterns in the sequence 
of A, B, C, D, A, B... in FIG. 7. The selection of 

10 the error diffusion technique is made by changing the 
weighting pattern to be used by the error diffusing 
section 52 in the predetermined sequence, thus 
minimizing moire and artifacts in the binary image, 
which did not appear in the original (multilevel) 

15 image . 

Moreover, as the error diffusing section 52 
proportionally distributes the occurred error to the 
plural unscanned pixels adjacent to the noteworthy 
pixel based on the weighting patterns as shown in (a) 

20 trough (d) of FIG. 7, the error diffusion technique 
changing section 53 changes the weighting pattern to 
another as selected at random from the above weighting 
patterns, thus minimizing moire and artifacts in the 
binary image, which did not appear in the original 

25 (multilevel) image. 

FIGS. 16 through 19 are graphs each showing the 
luminance (£*) of color chips of the image halftoned 
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by the halftoning apparatus of the present embodiment; 
the luminance, which was measured for every 5% 
throughout the range of from 0% to 100%, is graphed 
in comparison with a conventional halftoning method 
5 in which the pixel value of the noteworthy pixel itself 
is binarized. FIG. 16 shows the relation between the 
wavelengths and the luminance for magenta; FIG. 17, 
for green; FIG. 18, for blue; FIG. 19, for composite 
black (equal in RGB levels) . In FIGS. 16 through 19, 

10 the solid line shows the result of the present 
embodiment; the broken line, the result of the 
conventional art. 

As shown in FIGS. 16 through 19, the luminance 
falls near color chip No. 11 (50%) in the conventional 

15 halftoning method, whereas no such drop of the 
luminance appears in the halftoning method and 
apparatus of the present embodiment, guaranteeing a 
smoother halftoned image. 

When the multilevel image to be binarized is 

20 constituted with a plurality of images (e.g., a color 
image composed of a plurality of color separations such 
as cyan, magenta, yellow, and black) to be halftoned, 
sharing a substantially identical profile, the 
pixel-on-prof ile detection section 54 may 

25 discriminate whether or not the noteworthy pixel is 
a pixel on the profile for only one of the plural images ; 
and the result of this discrimination is temporarily 



stored in the RAM 16 or the magnetic disc 26, and the 
profile information (location or direction of the 
profile) may be used in halftoning all the images. 

All the plural images can be thereby halftoned 
without the need for the pixel-on-prof ile detection 
section 54 to perform the pixel-on-prof ile detection 
for the remaining images, thus speeding up the 
halftoning process. 

The plural images constituting the color image 
should by no means be limited to cyan, magenta, yellow, 
andblack as described above, but they may be red, green, 
and blue; other color separations also are applicable. 
Also any combination of colors constituting a color 
image is applicable. 

At that time, the binarizing all the images may 
be performed after completion of the profile detection 
or simultaneously with it. 

Further, in the halftoning apparatus of the 
present embodiment, as the multilevel image is scanned, 
the pixel-on-prof ile detection section 54 calculates 
the profile value E to discriminate whether or not the 
noteworthy pixel is a pixel on the profile, 
simultaneously with the binarizing performed by the 
binarizing section 51, thus reducing the time needed 
for halftoning. 

Since the halftoning method of the present 
embodiment is suitable for an parallel processing 
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operation, parallel processing using plural CPUs can 
be realized, thus speeding up the halftoning processes 
of a huge amount of image data . 

When the pixel-on-prof ile detection section 54 
5 carries out the profile detection, any one technique 
may be employed as selected from (1) using a Laplacian 
filter, (2) using a Prewitt filter, and (c) making 
addition and subtraction on the pixel values of the 
noteworthy pixel and the adjacent pixels. 
10 Furthermore, the present invention should by no 

means be limited to these embodiments, and various 
changes and modi ficat ions may be made without departing 
from the gist of the invention. 

For example, in profile detection by the 
15 pixel-on-prof ile detection section 54, other 

techniques than those used in the present embodiment 
may be used with the same results. 

And to calculate the profile value E by making 
addition and subtraction on the pixel values of the 
20 noteworthy pixel and the adjacent pixels, X ± is 

multiplied by 256/768 (=1/3) . The number by which X ± 
is to be multiplied should by no means be limited to 
this illustrated example "256/768 (=1/3)". 

Further, in the above description, the pixel 
25 value of the noteworthy pixel A is calculated directly 
from the multilevel values of pixels other than the 
noteworthy pixel A. This invention, however, should 
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by no means be limited to this, and the multilevel values 
of pixels other than the noteworthy pixel A are 
respectively converted into binary values before a 
two-dimensional digital filter is applied to errors 
accompanying the conversion. The resultantly 
obtained value is estimated to be the error of the 
noteworthy pixel that are to be diffused. 

The present invention should by no means be 
limited to the above-illustrated embodiment, and 
various changes or modifications may be suggested 
without departing from the gist of the invention. 
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